Central limit theorem for an additive functional of the 

fractional Brownian motion 



Yaozhong Hu* David Nualart^ and Fangjun Xu''- 
Department of Mathematics 
University of Kansas 
Lawrence, Kansas, 66045 USA 



Abstract 

We prove a central limit theorem for an additive functional of the d-dimensional frac- 
tional Brownian motion with Hurst index H G (j^, using the method of moments, 
extending the result by Papanicolaou, Stroock and Varadhan in the case of the stan- 
dard Brownian motion. 
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1 Introduction 

Let {^B{t) = {B^(t), . . . , B'^(t)),t > 0} be a (i- dimensional fractional Brownian motion (fBm) 
with Hurst index H G (0, 1). If Hd < 1, then the local time of B exists (see, for instance, 
[4, 5, 6]) and can be defined as 



Lt{x) = [ 6{B{s) - x)ds, t>0,xe 
Jo 



where S is the Dirac delta function. The above local time is jointly continuous with respect 
to t and X (see [4]). For any integrable function / : M'' — )■ M one can easily show the following 
convergence in law in the space C([0, 00)), as n tends to infinity 

(n'"'-^ j^' f{B{s))ds,t>Q) 4 (^Lt{0) J^J{x)dx,t>Oy (1.1) 
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In fact, making the change of variable s = nu, and using the scahng property of the fBm 
we see that the process (^n^'^~^ j^^ f{B{s))ds , t > j has the same law as 

rt r r 



^ / f{n''B{u))du = n'"' f {n" x) Lt{x) dx = / f{x)Lt{n-''x)dx, t>0. 

Jo jRd 

From here it is straightforward to verify (1.1). 

If we assume that J^af{x)dx = 0, then we see n^'^~^ J^* f{B{s)) ds converges to 0. It 
is interesting to know if there is a /3 > Hd-1 such that J^' f{B{s)) ds converges to 
a nonzero process. This will be proved to be true. In order to formulate this result we 
introduce the following space of functions. Fix a number /3 > and denote 

= [f eL^{W^): I \f{x)\\xfdx <oo and / f{x)dx = Q\. 

For any / G Hq^ by Lemma 4.1, the quantity 

f{x)f{y)\x-yfdxdy (1.2) 



is finite and nonnegative. The next theorem is the main result of this paper. 
Theorem 1.1 Suppose < H < ^ and f G . Then 

(n^^ j^' f{B{s))ds,t>Q) A (v/C^||/||^_,mLt(0)),t>0 

in the space C([0, oo)), as n tends to infinity, where " — ?■ " denotes the convergence in law, 
W is a real-valued standard Brownian motion independent of B and 

CH,d = -—J / w """{l- exp(-— ^)) dw = — -^r 



(27r)f7o ' 2w;2^^^^ (1 - iJrf)7rf V 2H 

Notice that -^^^ > Hd — 1 since H < When d = 1 and if = |, the above theorem is 
obtained by Papanicolaou, Stroock and Varadhan in [11] with Ci = 2. On the other hand, 

the constant CH,d is finite for any H > We conjecture that our result also holds for 
^TTTT < H < -,. But we have not been able to show our result in the case H < -i^. The main 

a+z d — a+1 

reason is that in the proof of Proposition 3.4 we need H > (see the Remark at the end 
of Section 3). 

In the critical case Hd = 1 the local time does not exist. For the Brownian motion case 
(if = I and d = 2), Kalhanpur and Robbins [7] proved that for any bounded and integrable 
function / : ^ M, 

1 z r 

/ f{Bs)ds A — / f{x)dx, 

log n Jo 27rV 

as n tends to infinity, where Z is a random variable with exponential distribution of param- 
eter 1. A functional version of this result was given by Kasahara and Kotani in [9], where 
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they also proved the second order results when J^2 f{x) dx = 0. The Kallianpur-Robbins 
law was extended to the fBm by Kono in [10], and the corresponding functional version was 
obtained by Kasahara and Kosugi in [8]. However, second order results for the fBm in the 
critical case Hd = 1 have not been yet proved. On the other hand, we refer to Biane [3] 
for some extensions of these results to the case of functionals of k independent Brownian 
motions. 

After some preliminaries in Section 2, Section 3 is devoted to the proof of Theorem 1.1, 
based on the method of moments. Throughout this paper, if not mentioned otherwise, the 
letter c, with or without a subscript, denotes a generic positive finite constant whose exact 
value is independent of n and may change from line to line. We use t to denote 



2 Preliminaries 

Let {-B(t) = (-B^(t), . . . , B'^(t)),t > O} be a (i-dimensional fractional Brownian motion with 
Hurst index H G (0,1), defined on some probability space {Q,J^,P). That is, the compo- 
nents of B are independent centered Gaussian processes with covariance 

E {B\t)B\s)) = ^(t^^ + -\t- 

The next lemma gives a formula for the moments of the increments of the process 
{W{Lt{0)) : t > 0} on disjoint intervals, where is a real- valued standard Brownian motion 
independent of B. 

Lemma 2.1 Fix a finite number of disjoint intervals {ai, bi] in [0, oo), where i = 1, . . . , N 
and bi < flj+i. Consider a multi-index m = (mi, . . . ^nij^), where rrii > 1 and 1 < i < N . 
Then 

N 

E(n[^(^^»(0))-^(^-(0))]") (2-1) 

1=1 

f n — ) f m- det{A{w))~'^^'^dw if all mi are even 

\t=l 2^ (2^)^ (mi/2)!/ Jn^iK.f'.]^ 

otherwise, 

where A{w) is the covariance matrix of the Gaussian random vector 

m ■ \ 

B(wl) : 1< i < N andl< k < —]. 



Proof. It is easy to see that when one of rrii is odd, then the expectation is 0. Suppose now 
that all rrii are even. Denote by J-"^ the cr-algebra generated by the fractional Brownian 
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motion B. Since ly is a standard Brownian motion independent of B, we have 

N 

El[[wiuxo))-wiL,M)r 

1=1 

N 

= E{E(n[m^..(0))-W^(La.(0))r|^^)} 

i=l 

N , N 

- [n5^]'^ni^.(o)-^.(o)r- 

= (n ."S^ .J / det(AH)-/2rf^. 

This completes the proof. ■ 

We shall use the following local nondeterminism property of the fractional Brownian 
motion (see [1]): for any = sq < Si < ■ ■ ■ < s„ < oo and Ui, . . . ,Un G M"', there exists a 
positive constant kn such that 

n n 

Var (^n, ■ {Bis,) - > A;^ J] - (2.2) 

4 = 1 1=1 

This can also be written as 



n n 



Var (^^Ui-5(si)) > A;^^^| (s^ - s^-i)'^. (2.3) 



i=l i=l j=i 

We claim that the law of the random vector [W{LbXO)) - W{LaXO)) : 1 < z < iV) is 
determined by the moments computed in Lemma 2.1. This is a consequence of the following 
estimates. Fix an even integer n = 2k, and set = {s e [0,t]'' : < si < S2 < ■ ■ ■ < Sk < 
t}. Let v4fc(s) be the covariance matrix of Gaussian random vector (i?(si), B{s2), ■ ■ ■ , B{sk))- 
Then the local nondeterminism property (2.2) implies that 

1 

(det Afc(s))~^ < c'' Y\_(^i ^ Si_i)~^'^ for some constant c. (2.4) 

i=l 

As a consequence of (2.1) and (2.4), 

E[WiL,{0))Y < c'nl [ sr^'^(s2-si)"^'^---(s,-s,_i)-^'^d5 

T{k{l- Hd) + 1) 

Therefore, E [iy(L((0))]" is bounded by c^'n!/r(/c(l — Hd) + l), and this easily implies the 
desired characterization of the law of the increments of the process {W{Lt{0)) : t > 0} on 
disjoint intervals by its moments. 
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3 Proof of Theorem 1.1 



By the scaling property of the fractional Brownian motion we see that, as processes indexed 
by t > 0, 

/ f{B{s))ds = / f{n"B{s))ds. 

Jo Jo 

Therefore, it suffices to show the theorem for the continuous process 



F4t) ■.= n-^ f f{n''B{s))ds. 
Jo 

The proof of Theorem 1.1 will be done in two steps. We first show tightness, and then 
establish the convergence of moments. Tightness will be deduced from the following result. 

Proposition 3.1 For any < a < b < t and any integer m > 1, 

E[(F„(6)-F„(a))'i <c((6-a)^-^'^ / \f{x)f{y)\\y\^-Uxdyy, 

where C is a constant depending only on H and m. 

Proof. Define 

D= [se M^"" : a < si < S2 < ■ ■ ■ < S2m. < h] . (3.1) 
Using the following identity for / G 

/(^) = 7^ / I e~^^-yf{y)dyd^, 

we then have 

2m 



(F„(6)-F„(a)) 

(2m)!n'"(^+^'^)E /" Y[ f{n" B{si)) ds 

1=1 



/„ \ , « « 2m „ / T 2m 2m 

[2m)l 



(27r)2' 



2md 



(2m)! 



(27r)2™'^ 

2m 

X 

i=l 



Imd 



D 



n /(^^) y^^^^ ( - 2^^^ ( E ■ ^(^^)) - ' E ^) 

2m „ ^ _ 2m \ 

n/(y.) / exp( --Var(5^e.-i?(s.)) ) 



n(exp(-.^)-l)rferfsdy, (3.2) 



where in the last equality we used the fact that Jj^^ /(x) c/x = 0. 
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By the local nondeterminism property (2.3), with the convention Sq = and ?72m+i = 0, 
we can write 



E 



2m 



(F„(6)-F„(a)) 

2m 

y ■ ■ 

X Y\ \ e~''^^ — l\d$^ds dy 
1=1 

f. 2m 



exp - 



2m 2m 
i=l j=i 



-cin 



m{l-Hd) 



md J 

2m 



i=l 



exp 



X 



2m 

I |2/ \2H 

-J- 2^ m {Si - 



dr] ds dy, 



where we made the change of variables rji = Y2'^=i for i = 1, . . . , 2m in the last equality. 



Let Xi = rii{si — and Ui = Si — for z = 1, . . . , 2m. Then 

"(F„(6)-F„(a))''" 

m{l-Hd) 



E 

<Cin 



2m 



2m 



n 1/(^^)1 (n 



U. 



-Hd 



exp 



2m 
i=l 



X 



2m 

TT exp ( ■ f ^^^^ — -tf) ) — 1 (ix (iw (iw . 



where X2m+i = and the integral on [a, 6] x [0,6 



12m-l 



means that Ui G [a, h] and 



G [0,6 — a] for 2 = 2,..., 2m. 

Let ^KuXx^ . . . , y^Kff^2m be independent copies of a d- dimensional standard normal 
random vector and ^2^+1 = 0. Then the above inequality can be rewritten as 



E 



2m 



< c2n-(^-^'^) / / ni/(^« 



(F„(6)-F„(a))' 

2 m ^ ^ 

xE ( TT I exp (i^^^-n—rr^ — L^^-^—^) — l| ) dudy 



j=l 



(3.3) 



Notice that 



11 < 2 A ( |e """^i - l| 



For each factor in the product inside the expectation in (3.3), we choose the upper bound 2 
when i is even and the upper bound 



e '+1 



l| + |e - l| 
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when i is odd. Thus, we have 



E 

<C3 



(F„(6) - F„(a)) 



2m. 



E 



R2md J[a,6]x[0,6-a]2'"-i 

H2i-l-^2i .H2i-l-^2i-l 



m ^ 

i=i ^ 



g " "2i-l 



Since the above random factors are independent, we have 



E 

<C3 



- F„(a)) 



2m 



R2m<i J[a,;,]x[0,6-a]2 



i=l 



X E ( |e """"2? 



^l/2i-l--^2i-l \ 

l| + |e """"^-1 - l| j J dudy. 

With the change of variables x = y2i^i, y = y2i, u = U2i and v = U2i-i, the above inequahty 
can be rewritten as 



E 

<C3 



(F„(6)-F„(a))^' 

n^~^'^\f{x)f{y)\{uvy^'^E(^\e ^^ - l| + \e'^ - l|) dudvdxdy 




b rb—a 



a Jo 



X 2 




fe— a pb—a 



m—1 



n^~^'^\f{x)f{y)\{uvy^'^E(^\e'^ dudvdxdy j . (3.4) 



'0 ^0 

Notice that J^^u'^^du < 04(6 — aY~^'^. Then, by Lemma 4.2, we obtain 

fb nb—a 




a Jo 



n^~"'^\f{x)f{y)\{uv)-"'^E (^\e'^^ - l| + - l|) dudvdxdy 

<c,{h-af-'"' f \f{x)f{y)\\y\^-''dxdy (3.5) 

JM.2d 



and 




b— a fb—a 



< C6(& — a 



^0 



^ |/(a;)/(2/)|(H E ( \e'^^ - 1\] dudvdxdy 



\fix)f{y)\\y\H '^dxdy. 



(3.6) 



Now Proposition 3.1 follows from (3.4), (3.5) and (3.6). 



Now we prove that the moments of Fn{t) converge to the corresponding moments of 
WiLtiO)). 
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Fix a finite number of disjoint intervals (aj,6j] witli i = 1,...,N and 6j < aj+i. Let 



m = (mi, . . . , mAf) be a fixed multi-index with nii E N for i = 1, . . . , N. Set YliLi — l^^l 
and YliLi ~ need to consider the following sequence of random variables 

N 
i=l 

and compute lim E(G„). Notice that the expectation of G„ can be written as 

n— >oo 

/N rrii 
l[l[f{n''B{si))ds), 

where 

Dm = {se : ai < s\ < ■ ■ ■ < si^^ < bi,l < i < N] . (3.7) 

Here and in the sequel we denote the coordinates of a point s E M'™' as s = (s*), where 
1 < i < N and 1 < j < rrii. 

For simplicity of notation, we define 

Jo = ■.l<^<N,l<J< nii}. 

For any (^l, j'l) and («2, ^2) £ </o; we define the following dictionary ordering 

(«l,Jl) < («2,J2) 

if Zi < Z2 or ii = 12 and ji < j2- For any in Jo? under the above ordering, {i,j) is the 
(X]fc=\ "^fc + i)"th element in Jq and we define = X]fc=i "^fc + J- 

Proposition 3.2 Suppose that at least one of the exponents nii is odd. Then 

hm E (Gn) = 0. 

n— J-oo 

Proof. The proof will be done in several steps. 

Step 1 Using similar argument as in (3.2), we obtain 

1 „ „ „ N rrii 

m! \in\(l-Hd) 



(27r) I""!"" Jui-i'* Jd^ ^ f^-^ fj[ 



X exp - 2Var (J:^:^; ■ 5(^})) - ^ E E 

WWm 



m! |m|(l-gd) 
7?7, 2 



i=l J=l i=l jr'=l 

N nii 



(27r)l'"l"' 

/Ui \ iv nii 



jTtj \ Af mj 



i=l j=l ' i=\ j=l 
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By the local nondeterminism property (2.3), with the convention Sg = s^.^^ for 2 < i < N 



where we used the fact J^j_ f{x) da; = in the last equality. 

h 

0, 

N rrii 

v-(EE?rB(^;))£«"5:5:| E a 

i=i j=i 



and Sq 



N nii 



i=l j=l {l,k)>{i,j) 



\2H 



Let F{y) = Hili Ilj^i /(^/j) and make the change of variables 77] = Y.{e,k)>{i,j)^k ^^i 
1 < i < N and 1 < j < rrii. Then we can estimate E[Gn) as follows: 



\m\{\-Hd) 

|E(G„)| <cin ^ 



N mi 

nn 



'|d J D 



exp 



|F(|/)|exp 



|m|(l-gd) 



„ „ N mi 

ij \ny)\{\m{s)-su 

jRl'^ldJDm \ = 1 ,=1 



TV mi 



(i?7 ds dy 



-Hd 



X 



AT mi 

i=i i=i 



— L 



dsdy, (3.8) 



where X^^^^ = 0, X^^+i = Xj+i for 1 < i < - 1, and ^EJ^Xj {1 < i < N, 1 < j < rui) 
are independent copies of a d- dimensional standard normal random vector. 



Denote the expectation in (3.8) by /. That is, 



/ 

where 



N mi 



nnn 

i=i j=i 



exp 



— L 



E 



n ^^^y 

(«j)eJo 



y)-x: 



for 1 < i < and 1 < j < rrii. 

Notice that the random variables lij for e Jo are dependent. We are going to 

choose a proper subset of Jo in the following way. Assume that rrii is the first odd exponent. 
Then we choose all the factors lij such that < ^{i,me) and is odd. Then, we 

choose all the factors lij such that > ^{i,me) + 1 and is even. Notice that 

all these factors are mutually independent and they are also independent of the product 
Ie,mile+i,i- The lack of independence of the two factors Ie,me and Ie+i,i will be compensated 
by the fact that the integral of {s{'^^ — s^^^)~^ is finite for any (3 < 2, because we have the 
constraint s^^^ < be < s{~^^. To make this argument more precise, let us define 

Je = Je,i U Je,2, 

where 

Je,i = G Jo : #(i,j) < #(^,m^) and #(i,j) odd} 
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and 



■h,2 = {{hj) e Jo : j) > + 1 and j) even}. 

Notice that lij < 2 for all G Jo- Then 



I < < 



iie = N. 



Step 2 We first consider the case i ^ N. In this case, the number of elements in J^ is 



^1 — 1 and 



\m\(l-Hd) 

|E(G„)| <C3n -2 



J Ml JDrn {i,j)eJo 
xE{le^raJe+i,i) n ^ikj)dsdy. 

In the last inequality, we used the fact that all random variables Ii^mili+i,i and /jj for 
(z,j) G J I are independent. 

Since je'^^^-^^) - l| < |e'^i - 1| + je'^^ _ i| for all ^i, 22 e M, 



Hi A. -J .^H 



I -"V-l) — 1|) 
X (^|e " ''l ) _ 1| 4- |e " 1=1 =m^J _ 



°1 

-1 



Notice that X2'^^,Xf'^^ and are independent. As a consequence, we can write 



where 



E < I[ + /. 



2) 



and 



/j^ = E |e '1 ""f' 



+ E le'"''^^™^-^'"^.-!^" - 1| E |e'"^(4+^-V^)« - 1| 
+ E|e * ™^ - 1|E 



(3.9) 



Therefore, 



|E(G„)|< cgn^^^^ / |F(2/)|(Gi + G2)Jy, 

|m|ti 



(3.10) 
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where 

Gk,n = j (y 



for A; = 1,2. 
We claim that 

In fact, making the change of variables = — s*„j^ for all (i, j) G Jo, and defining 



l^i^i+J-^'^E(Ji,,) if G J, and ^ (iV,m;v); 



Af \-Hd 



E(/Ar,mjv) if («, j) e and = {N,mN), 



and = (f2'^^'yi'^^t'L) ^'^ ^i, we obtain 



Here we used Lemma 4.2 in the last inequality f-^] + 1 times. 

For any f3 G [0, 1], we have |e^^ — 1| < C/3 for all 2; G M. Recall the definition of in 
(3.9). We then have 

<crn-'''%Lf\yi^risi^'-sLr''''- 

So 



Define J^,3 = G Jo : j) < #(£,mf)} and 

= {«i < 4 < ■ ■ ■ < 4, < 1 < ^ < ^ - 1; < ^1 < ■ ■ ■ < 4,-1 < ^4- 
Integrating the above integral with respect to for (i, j) G J^,2 and using Lemma 4.2, 
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where #J^,2 is the cardinahty of J and 



-Lt-i ■^"■e+i 

We observe that, if 1 - Hd < 2Hf3 < 2 - 2Hd, 



I 



nil/ ^ "iTT-t m£ — i/ "'"mi 



As a consequence, 



^ n - n 

= ci3n-^^^n([^l-i)(^^-i)bij/^b^r n l^^-l^"'- (3-12) 

Choosing (3 = jj — d in (3.12), 

G,,„<ci3n([^]+^)(^^-^)|ylj^"1y^^|^-^ n l^^-l^"'- (3-13) 

Substituting (3.11) and (3.13) into (3.10) yields, 

|E(G„)|<Ci4ri-^ / iFivW^^l^-'lyi-^'l^-" J] W^'^dy. (3.14) 



Step 3 Now we consider the case £ = N. In this case, Je = Ji^i and 
\E{G^)\<c,,n''^ f I \F{y)\ \[{s)-s)_. 



Hd 

(i,i)eJo 

xE(/;v,m^^) n ^i.ki)dsdy. 



(*,i)6^jv.i 

Define Jat^s = {(i, j) G Jq : #(«, j) < #(A^,"^7v)} and 

= {oi < si < ■ ■ ■ < s^^ < 6i, 1 < i < - 1; a£ < sf < ■ ■ ■ < < &jv}. 
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Integrating the above integral with respect to and using Lemma 4.2 yield 



|E(G„.)|<c,en """-T-' / / \F{y)\\y. 



TTtArl 



X n I 

(lj)eJiV,3 



Using similar arguments as in Step 2, 



\E{Gn)\<c^^n- 



|m|d 



(j,i)eJjv,i 



(3.15) 



Step 4 Recall that / G H^" . Then, from (3.14) and (3.15), we see that |E(G„)| is 
bounded by a multiple of n ~ . Our result now follows from taking the limit. ■ 

In the sequel, we consider the convergence of moments when all exponents are even. 
Recall the definition of in (3.7). For 1 < £ < and 1 < < ^, we define 



^2k ^2k-2 



< ^2k ^2k- 



The following result tells us that the integrals over the domain Ol do not contribute to the 
limit of the moments. This result will play a fundamental role in computing the limits of 
even moments. 

Proposition 3.3 For any 1 < £ < N and 1 < k < ^, 



„ N rrii 

Jim n^^^^E ( / ^ n n f^^^'Bis))) ds) = 0. 



Proof. Using the arguments and notation in the proof of Proposition 3.2, we obtain 

„ A'' nrii 

^{l\{\{f^^''B{s)))ds] 



\m\(l + Hd) 

n 2 



fe 1=1 j=i 



„ A'' rrii 

Kcn'^i / iFfe)i(nn(»;-»;-i: 



X 



TV mi 

i=l j=l 



ds dy, 



where X^^^^ = 0, X'^^^^ = X|+^ for 1 < z < iV - 1, and ^^-^j {I < i < N, 1 < j < rrii) 
are independent copies of a dimensional standard normal random vector. 

We make the change of variables = — for all G Jo- The integral domain 
Ol becomes 

Di={ve R^f : ai < J2 < ^2.-1 <<k}- 

(ij)GJo 
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For G Jo, define 



Then 



|m|(l+Hd) 

n 2 



N mi 



< c\ n~ 



f [ F{y)Ef H {v])~"'h^dvdy. (3.16) 



|m|(l-fl"d) 



Next we estimate the expectation in (3.16). We are going to use an argument similar to the 
one used in the proof of Proposition 3.2, based on the selection of some factors in the above 
product. Here, the dependent product that will play a basic role will be l£2k h2k-i, due to 



the definition of the set O^. Define 



f - f IJ f 



where 



4.1 = {(hj) e Jo : #(z,j) < #(^,2A;-2),#(z,j) odd}, 

4.2 = {(hj) e Jo : #(2,j) > #(£,2A:),#(^,j) even}. 



Since all exponents rrii are even, the number of elements in J^ is 



|m|-2 _ |m| 



1. From 



the definition of Jjj, we know that random variables Ii,2kle,2k-i and lij for {i,j) G J^ are 
independent. Then 



\m\{l + Hd) 

n 2 



N mi 



< C2n' 



I I F{y)E{I,,2kh2k-i) n H)'''' n ^ik,)d^'dy. 



|m|(l-gd) 



(«j)G^O 



For G Jf and ^ {N,mi\f), define 



and Oat, 



[0,biv]2 



2, we 

a,,,dv]dv}^,<C2n''''-'\y]\i^' 



N ^ From Lemma 4.2, we obtain 



for all G Jf. Therefore, 



|m|(l+Wd) 

n 2 



N rrii 



fe i=l j=l 



< C2n 



l-Hd 



F{y) n \y}\^~'4dy, (3.17) 



(*j)6Jf. 
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where 



Jo Jvi,_, Jo 



Notice that |e'*^^^ _ i| < jg'zi _ i| + |e^^2 _ ]^| f^j. gj[ Z2 E M. Using the independence 
of X2k^^,X2k and X^^j^i, we obtain 



e. 

fc,25 



E (/£,2fc-i/£,2fc) < 4,1 + ^: 

where 

4 1 = E \e^^^^^^ -l\E\e ""^^^k+i)" _ i| + e |e - 1| E |e ^^^^^^^ - 1| 



y2k-i'^2k-i y2k'^2k 
+ E|e ""("^.-1)" _i|E|e'""("^.)" -1| 

and 

4,2 = E(|e^""(''^.)" _l||e^"«(^^.)« _1|). 

Now we have 

4' = 4' 1 + ^2, (3.18) 

where 

= III ('^2A:+l'^2fc'^2fe-l) ^k.idV2k+ldV2kdV2k-l 
Jo Jvi,^, Jo 

for i = 1,2. By Lemma 4.2, 

/li<C3n-(^-^^)bL-#-1l/L|--'^. (3.19) 
For any /3 G [0, 1], we have |e'^ - 1| < C/sl^;]'' for all 2 e M. Then, 

Therefore, ii 1 - Hd < 2Hp < 2 - 2i7rf, 



Jo 



2k-l 





v\k^,r"\AK)"'''~''''{A,^ir^'^ dv^2k^, dvi, dvi,_, 

6jv 



Choose (5 = ^ii^ 



^ ^ 3(1--H"d) , ^ , 3(l-gd) , ^ , 3(l-gd) 



IL < cjn ^\yk-i\^^\yk\^^- (3.20) 
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Substituting (3.19) and (3.20) into (3.18), we obtain 



„ 3(1-Hd) . „ 3(l-fl"d) , „ 3(1-Hd) , „ , 1 J, « , 1 

4 < csn — ^{\yL^^\^^\yk\^^ + \yik^i\^~'\yk\^~')- (3.2i) 

Our result now follows easily from (3.17), (3.21) and the assumption / G Hff ■ 

Before proceeding with the proof of the convergence of moments when all exponents rrii 
are even, we would like to make some heuristic remarks in the simple case E(F„(t)^), that 
might help the reader to understand our technical computations. We can write 

E(F„(t)2) = 2^^+^^^ r Te (/(n^5(si))/(n^i?(52))) ds.ds,. 
Jo Jo 

Making the change of variables ui = n{s2 — si) and U2 = S2 yields 

E (F„(t)2) = 2n^'^ E (^f{n"Biu2 - ^))/(n^i?(n2))) dm du2 

= 2n^'^ 11 I f{n^x + y)f{n^x)pn{x,y)dxdyduidu2 

Jo Jo JR^d 



t pnu2 

Jo 



'^11 I f{x + y)f{x)pn{^,y)dxdyduidu2, (3.22) 



where Pn is the density of the 2(i-dimensional random vector {B {U2) , [B {u2 — ^) — 
B{u2))). Then, as n tends to infinity, u^^n^ {B{u2 — ^) — B{u2)) converges in law to 
a (i-dimensional standard normal random vector independent of B. So, formally we obtain 
that the expectation E (F„(t)^) converges to 

e( [ 6{B{u2))du2) [ [ f{x + y)f{x)u^^'^(e~^ -l)dydxdui, 

^Jo ^ Jo ilR^d V / 



(2vr: 



which is equal to C*H,d||/|| j__^E (H^(Lt(0))^) . Notice also that we have been able to add the 

term —1 because the integral of / is zero. In conclusion, the term B{u2) appearing in (3.22) 
contributes to the local time at zero whereas (^B{u2 — ^) — B{u2)) becomes independent 
of B in the limit and its expectation contributes to the constant Cha- When computing an 
even moment, for each couple of consecutive factors we will observe this phenomenon. 

The main difficulty to make this argument rigorous is to compute the limit of the integral 
of the density p^ over the interval [0, nu2\- To overcome this difficulty we will first integrate 
on a compact [0,K], and then show that the integral over [K,nu2\ converges to zero as K 
tends to infinity, uniformly in n. However, this convergence holds only if we integrate over 
[K^ n|2j^ which is fine because Proposition 3.3 implies that the integral over [^^,nu2\ tends 
to zero as n tends to infinity. 

Consider now the convergence of moments when all exponents are even. On each 
portion of the coordinates ai < s\ < ■ ■ ■ < < hi we make the following change of 
variables: 

U2k = -slk and M2fc-i = nis2k - S2fc-i)> where 1 < A; < mi/2 
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with the convention Uq = Sq = a.j. The idea is to couple each variable with an odd subindex 
with the next one. In this way we obtain 



N -t 



E = m!n^E ( / H H ][n'' B{^^^))j{n'' B[v^^^ - ^)) dv) , (3.23) 



m i=l k = l 

where K and -D" cirG cis follows: 

K 

and 



\m\Hd 



0-^ = {ue : a, < ^ <<<■■■<<, < 

< ?4;._i < n(u*2fc - W2fc-2), 1 < ^ < 

We compute the expectation (3.23) in the following way. Define the m- dimensional 
Gaussian random vector X{u) by 

Xi,{u) = Biul,) and X^,„^(«) = n^(i?(«^, - ^) - 5(«y), 

where 1 < k < The covariance matrix and the probability density function of the 
Gaussian random vector X{u) are denoted by Qn{u) and 

|m|d 1 , 1 rp^ 

Pn{x) = (27r) 2~(det(5„(n))~2 exp [xQn{u)~ x ) , 
respectively. With the above notation we can write 

r r ^ ^ 

E (Gn) = m!n^ / / H H f (^"^ ^2k) f {n"" x\^ + x\^^,)pn{x) du dx. 

Making the change of variables yj = n^x'j if j is even, and = Xj if j is odd, we then 
obtain 

N 

E (G„) = m! / / n n /(2/k)/(2/2. + yk-i) Pn(y(^)) c/n c/y, (3.24) 
where ?/](n) = n~^yj if j is even and i/](n) = if j is odd. 

Proposition 3.4 Suppose that all exponents rrii are even. Then 

I ^ 

hm E(G„) = CZ WfWilMU^WiUXO)) - W{L^XO))r). (3.25) 



i=l 



Proof. The proof will be done in several steps. 

Step 1 Notice that we can find a sequence of functions /at, which are infinitely differen- 
tiable with compact support, such that J^^ In^x) dx = and 



lim 

N^oo 



\f{x)-fN{x)\ {\x\Tr-^Wl')dx = 0. 
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So, by Proposition 3.1, we can assume that / is infinitely differentiable with compact support 
and Jjgd f{x)dx = 0. 

The equation (3.24) can be written as 



E = m! / / F{y) pMn)) du dy, 



(3.26) 



where 



N -T- 



^(^)=nn/(^2j/(y2fc+2/k-i)- 

j=i fc=i 

Let us compute the hmit of the density Pn{y{n)) as n tends to infinity. We spht the ran- 
dom vector X{u) into two random vectors X{u) = {Y{u), Zn{u)), where Y{u) contains the 
components of X{u) with even sub indices, and Zn{u) contains the components with odd 
subindices. That is, Y{u) is an ^^^-dimensional random vector, such that Yf^{u) = B(u2f,) 
for 1 < z < and 1 < k < We denote by A{u) the covariance matrix of Y{u), which 
does not depend on n. On the other hand, the covariance matrix between the components 
of Zn{u) and Y{u) converges to the zero matrix, and the covariance matrix of the random 
vector Z{u) converges to a diagonal matrix with entries equal to (M2fc-i)^^; 1 < ^ < A^, and 
l<k<^. Therefore, 



N 

hm p„(y(n)) = (27r)-^(det A(«))-^ nn("2fc-i)"'''exp 

n^oo J. J. J. J. 

1=1 fc=l 

On the other hand, the region D"^ converges, as n tends to infinity, to 



\V2k-\\ 

2u 



2k-l 



|m e mI™I ■.ai<ui<---<ui^^<bi;0< U2k-i < oo; 1 < A; < 1 < i < A^j. 
Notice that we can add a term —1 because /j^^ -F(y) dy\j^_^ = for any i, k, and 

du = -|y^fe_il^""' 



-Hd 



e 2u 



-Hd 



1 — e 2n 



du. 



Therefore, provided that we can interchange the limit with the integrals in the expression 
(3.26), we obtain 



where 



hm E(G„) = m!2-V(2^)-^C^^J|/||fi 



|m| 



{detA{w))-^dw, (3.27) 



Or 



Oss. = <w eR 2 : tti <wl < ■ ■■ < w\ <hi,l<i< N\, 



and A{w) = A{u) with the change of variable wl = u^i^- Finally, the right-hand side of 
(3.27) can also be written as 



N 

n 



-CI 



{27ry— {detA{w)y 2 dw, 



(3.28) 



and, taking into account Lemma 2.1, this would finish the proof. 

Step 2 In order to justify the passage of the hmit inside the integrals, we decompose the 
region into two components as follows. For Ji" > 0, we define 



2k "'2k-2)'i ^ < 



and Dl^^^, = n^- D^^^,. Then, E (G„) = I^.k + IIk. where 



m! 



F{y)Pn{y{n))du dy, 
F{y)Pn{y{n))du dy. 



The region ki '^^ uniformly bounded in n and we can then interchange the limit and the 
integral with respect to m, provided that we have a uniform integrability condition. To do 
this we need the following estimate of the density Pn{y{f^))- 

For any ^ G M'™' with components (^j), l<i<iV, l<j< mj, we can write 



i=l ^ 7 = 1 k=l 



u 



u 



2k 



1=1 ^ j=i 

N -t 

i=l fc=l (e,2j)>{i,2k) 
N -w- 



2fc-l ' 

n 



2fc-2. 



2fc 



2k-\ ■ 



i?(4.-2)) 



+ E E ( E 4- - ^""^2.-1 ) ■ ( B{v},,) - B{u 

i=l k=l {e,2j)>{i,2k) 



u 



2k 



2fc-l ■ 

n 



Here we have used the ordering (£, 2j) > (z, 2A;) if £ > -i or £ = z and j > k. 
By the local nondeterminism property (2.3), 



Var(e,X)>A;H 



N -t 



EE I E «^,r-K " 

i=l fc=l (e,2j)>{i,2k) 



2k~l i \2i? 



N 



+EEI E fi (^) 

i=l k=l {e,2j)>{i,2k) 



2 2H 



AT 



j=l fc=l 



ha/fcl (U2k u 



n 



2k~2 



2H 



i=l k=l 



(3.29) 



19 



where we have made the change of variables 

Vik = iij and r/^fc_i = ^^^^.i- 

(£,2i)>(i,2fc) 

This imphes that 

(detgO"^ = (2vr)-^ / exp(-ivar(e,X))(ie 



I ml d 



< (27r)-^ / exp ( - ^Rij]) ) dr] 



N 



i=\ k=l 



n 



Therefore, 



Pn(z/(n)) < C2 n n K.-i)^'''(^2. - ^ - ^'2.-2)' • (3.30) 

i=l k=l 

As a consequence of (3.30) and the inequahty (4.5) in Lemma 4.5, 
/ Pn{y{n))du 

r ^ ~ ^ —Hd 

<csi nnK-i)"^'(^.-^-«2.-2)" du 



i=l fc=l 



< C4, 



where C4 is a constant independent of n and y. Thus, taking into account that the function 
F{%j) is integrable, by the dominated convergence theorem we obtain 



n— >oo 



i! / F{y)( hm / pn{y{n)) du)dy. 

m, A ,1 

On the other hand, again by (3.30) and Lemma 4.5, there exists p > 1 such that 

sup / \'Pn{y{n))\^ du < 00^ (3.31) 

m,X,l 

which imphes 

lim/ii^ = m!/ / F{y) lira ^ ^{u)pniyin)) dudy. 
With the same notation as above we get 

hm/^^ = m!(27r)-^f / (det v4(tf;))-^rfy;') (3.32) 

"2" 



X 



/ n n f /(2/2.)/(2/2. + 2/2.-1) / n-^'^(e-^ - 1) du) dy. 
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The right-hand side of the above equahty converges to the term in (3.28) as K tends to 
infinity. 



Step 3 Now it suffices to show that 

n— ^oo 

First we observe that 



hm hmsupJ^^ = 0. (3.33) 



where 

So we only need to show that for all (i, k) 



m,K,i.k 



lim limsup / / F{y) pn{y{n)) du dy = 0. (3.34) 



As a consequence of Proposition 3.3, we can replace -DmXifc (3.34) with 

Dm,K,i,k = [ueD^; K < U2fc_i < ^ 1 

and just show that for all (i, k) 

lim limsup / / F(y) Pn{y{n)) du dy = 0. (3.35) 

K^oo n^oD JMlm|d Jd"'\^ . , 

To do this we need more refined estimates of the density Pn{y{n))- By Fourier analysis 

pMn)) = (2vr)-l-l'^ / exp f - ^Var (e, X> - . ^ f + ^2,^, ■ y^,^^ d^. 

j]Rim|d \ / .^^ \ n / 

We choose a set J of indexes of the form {i,2j — 1), where 1 < z < and 1 < j < 
For each index in J we introduce the operator 

A,2,„iF(y^^.„i) = Fiyl^_,) - F(0) 

and set Aj = Y[[i2j-i)ej Taking into account that the integral on the variable 

ylj-i is zero, we can replace Pniyin)) in (3.35) by Ajp„(?/(n)). Using (3.29), we obtain the 
following estimate 

/^jPn{y{n)) < C5 / exp ( - ^ V V [|r/^/(M2i " " ^2^-2)^^ 

jR\rn\d \ z .^^ .^^ L n 

i \ N ^ 

+ \rih - n"r,^2,-^\^r] ) JI 11 1^""^^-^"'^-^ " ^1 

^ ' i=ij=i 

JlR|m|d \ Z L 



cgn 2 / exp ( - — > J> J |r/2j.| (^Waj W2i-2 



lit ■ 



+ l^2,-i|(^) JJlUlle -l^r/. (3.36) 
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This shows that 

F{y)Pn{y{n)) dudy\ 



m.K.i.k 



^ / i=ij=i 

To estimate the right-hand side of (3.37), we first consider the integral in the variables 
u = ^2^, V = w = ?72fc and z = ?72/fc_i- Set ^2^-2 = % 2/2fc-i — V- That is, we have 

the integral 



uo J K 



We see that 
where 



x(|e '■"^ — l| + |e '^'^ 

/i(mo, y) = hiuo, bi, y) + /3(mo, y), 



I I 
X e~^^ — 1 (iw; (iu 



and 



, _ n(M-tio) 
•^2 r :^ 



h{uo,bi,y) = n 



-Hd 




uo J K 



X e"'^ — 1 dw dz dv du . 
Integrating with respect to w and using the inequality (4.4) in the appendix lead us to 



{u-uq y"'^ I e"^^^*^"^ \e~''^ -l\dzdvdu 



UQ J K 

Hd 



(n - uoY^^du 

< cK'-'"'-''^{h-uof''"'\yf. (3.38) 
In a similar way but with the application of (4.3) instead of (4.4) we obtain 



n(u-~U(\ ) 

v-^ e \e-'^ - l\ dw dv du 

J un 



< cni-^'^-^^(6, - u,f-^'''-''^\yf . (3.39) 
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Combining (3.38) and (3.39) gives 

h{uo, k, y) < cKi-^"-^/^(6, - a,)i-^'^|y|^ + cni-^'^-^^(6, - a,)2-2^^'^-^/^|^|^ . (3.40) 

Once this is done, we proceed to consider the integrals in the variables u = u^i, v = U2i_i, 

z = r]2i and w = ?72;_i with indices (j, /) 7^ (z, k). Set also ti2/_2 — ^0 and y2«-i — V- That is, 
we have the integral 



JuQ Jo ^]R2d \ Z \ n n // 

|e~';r7r - l| + |e" ^ - l| j dwdzdvdu. 



We can decompose this integral into two components: 

h{uo, bj, y) = h{uQ, bj, y) + h{uQ, bj, y), 

where 



X e~ ^ — 1 dw dz dv du 



[u — Uq j / e 2 V™/ |e ^ — l\dzdvdu 



and 



= „-™/7"'"""7 exp(-^|i(|.P(;i)- + |uf(«-.„-^)-)) 



'uo Jo 

X e~ ^ — 1 dw; (if 



y~Hd ^-^i^[u-uo~-) \e~'^ -l\dwdvdu. 
The inequalities (4.2) and (4.1) imply that 

h{uo,bj,y)<c [ \u - uo^^lyli^-" du 



-Hd\„,\^-d 

c(6,-«o)^-^1z/P 



\l-Hd\„.\jj-d 



<c{b,-a,f-^'\y\^-'^. 

The remaining integrals can be dealt with in the same way as l4{uo,bj,y). Thus the state- 
ment (3.35) follows. The proof is completed. ■ 



Proof of Theorem 1.1. This follows from Lemma 2.1, Propositions 3.1, 3.2 and 3.4 by 

the method of moments. 
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Remark Although the constant Cn^d is finite for H > -j^, the proof of Theorem 1.1 only 
works for H > The reason for this is that for any ?/ G M*^ 

(-f - ^"") ^« - (I)' - (-2!!!^)) • 

which is bounded by c\y\'^u~^^'^~^'^\ while, on the other hand, 

f exp f-^lel V^) \e'^-y -l\d^< c|y|w-^('^+i). 

So, any type of estimation procedure, like the one based on the local nondeterminism prop- 
erty used in this paper, will lead to an upper bound of the form u~^^'^~^^\ 

4 Appendix 

Here we give some lemmas which are necessary in the proof of Theorem 1.1. 
Lemma 4.1 Let < (3 < 1. If f E Hq , then ||/||^ given in (1.2) is well-defined and 

ml = [ \Tm)m-^-'d^>o, 

where J^f{^) denotes the Fourier transform of f and 



[ {l-cos{x-Om-^-'d^>0 



is independent x if \x\ 



Proof. For any x E S'^ ^ and any dx d orthogonal matrix Q, the change of variable ^ = Qr] 
yields 

I {I - co&{x ■ d^ = I (1 -cos((Q'^x) ■r/))|r/|-^-^dr/ > 0. 

This shows that (l — cos(x ■ Cj)\i\~^~'^d^ depends only on The substitution ^ = -j^?] 

gives us Jjgd (1 — cos(x ■ 0) 1^1"'^"'^ '^^ = Ci3^d\x\^ ■ Then an elementary result from Fourier 
analysis [12] yields 

\l = c~^' [ f{x)f{y)( [ {e^^^"y>^-l)\^r^~''dAdxdy 

= c-^' [ \j'f{om~^-'d^>o- 



Lemma 4.2 Assume that Hd < 1. Let X be a d-dimensional centered normal random 
vector with variance cr^. Then, for any n E N and y G M'^, there exists a constant c 
depending only on H, d and a such that 

^Hd~l\„.\^-d 



f .-E|exp(,|^)-l 



du < cn'"'~'\y\H' 
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Proof. It suffices to show the above inequahty when y ^ 0. Making the change of variable 
V = \y\ unu gives 



du 



n 



Hd-l\ 



dv 



< n 



Hd~l\ 



v^"'^(2Air^E\X\j dv 



cn 



^-Hd-Hp 



Lemma 4.3 Assume 1 - H < Hd < 1 and < P < 1 with HP < 1 - Hd. Then, for all 
n G N, there exists a constant c independent of n such that 

Jo Jm-i 

(4.1) 

/ {s--y du e-\^\'^''^^'"\l-e'^\di<cn-^h-"'^-''^\y\^-\ 

Jo ^ jRd. 

(4.2) 

where k is a positive constant. 

Proof. We ffist note that (4.1) follows easily from (4.2) by the change of variable. So, it 
suffices to prove (4.2). Denote the left hand side of (4.2) by /. We then have 

I = n~^'~^^ r is du [ e-\^\'&"\l-e'i^\d^ 











n 



where we made the change of variable C,{^)^ = x. 
By Lemma 4.2, we obtain 







n 

f oo 



f^_^yHd^H^^^Hd / e-^\'\'\l-e'T^\dxdu 



poo /• 

< c,s-^'^-"(^ / u-^'^ / e-"l"l'|l - e'T^ldxdu 
Jo Jr'^ 
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On the other hand, 



J ns n 

where f3i can be any constant in [0, 1] and we used H/3 < 1 — Hd in the last inequahty. Our 
result follows by choosing f3i = — d. This is possible because 1 — H < Hd < 1. ■ 

Lemma 4.4 For n E N, we assume that 1 — H < Hd < 1 and < K < Then there 
exists a constant c independent of n and K such that 

ns 

J K Jw 

ns 

n-""' r [s du [ e-«l«l'(^)'"|l-e^&|cie<cs-^'^A'i-^'^-^^|y|^ (4.4) 

where n and /3 are positive constants with 1 — Hd < HP < H A {2 ~ 2Hd). 

Proof. The inequality (4.3) follows easily from the proof of Lemma 4.3. We shall show (4.4). 
Denote the left hand side of (4.4) by /. Then we have 



I = n 



J K/n jRd. 

Jk 



^''\s-u)-"''u-"'-''Pdu. 

I K/n 



Since "f- < f and 1 — Hd < H(3, we have 

r^\s - uY'^'u-'''-''^ du < OS-'"' f'\-'"'-''^ du < cs-'"'{-y-"''-"^. 

JK/n JK/n 

Therefore, 

This proves the lemma. ■ 

Lemma 4.5 For any K > and n G N, there exist constants ci and 02 independent of n 
such that 

KAnu 

y-Hd^^ _ ly^'^dv < crK'-'^'u-''', (4.5) 
^ 
and 

fU 



V 



{u-vY"''dv = C2U^-^'"'. (4.6) 
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Proof. The inequality (4.6) follows easily from 



1 



Jo 
We only need to show (4.5). Notice that 

pKAnu p^Al 

Jo n' Jo 

If nu < 2K, then 

Jo 

If nu > 2K, then 
Therefore, 

pKAnu 
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